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1. Introduction 


Considerable interest has been attached to the research and develop- 
ment of free radicals and electronically excited metastable species for 
use in rocket propellants, because markedly superior performance is pre- 
dicted theoretically for then. It is well-known that the maxinun for 
specific impulse with the best conventional chemical propellants is in 
the neighborhood of 500 sec. On the other hand, values for the theoreti- 
cal specific impulse in the neighborhood and even well above 750 sec are 
predicted if it is possible to employ the ek 't promising free radicals or 
electronically excited metastable species (i.e., "metastables") . Such a 
50% Increase in specific impulse would engender a dramatic concomitant 
increase in payload and a decrease in the number of stages in space vehi- 
cles (Carter^®, Windsor^^) . 

This report provides an up-to-date survey study of free radicals and 
electronically excited raetastable species as high energy propellants for 

rocket engines*. Included as free radicals in our study are nascent or 
atomic forms of diatomic gases in addition to the highly reactive dia- 
tomic and triatomic molecules that possess one or more unpaired electrons. 
For use as a propellant, a free radical must be a stable (ground-state) 
quantum mechanical structure and have a rather high energy of reaction 
per unit mass, which implies that the free radical must be of relatively 
low molecular weight. Electronically excited metastables of possible 
Interest as propellants must have sufficient resistance to spontaneous 
electromagnetic decay (i.e., a sufficiently long radiative lifetime) in 
addition to having a high specific energy content. 


^Metallic hydrogen and vibrationally excited molecular species, also of 
interest as possible high energy propellants, are not discussed in this 
report. A recent description of current experimental research on the pro- 
duction of metallic hydrogen has been given by LubHn^. Noteworthy ex- 
perimental and theoretical results on vibrationally excited hydrogen, oxygen, 
and other molecules have been obtained recently and reported in Refs. 37, 

39, 54, 55, 76, 103, 109, 123, 160, 182, 195 and 207. 
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Free radicals and metastables that nay offer possibilities as pro- 
pellants are discussed in Sections 2 and 3. In section 4 it is shown 
that propellants composed of any of the free radicals H, N, 0, CH, BH, 

NH, BH 2> CH 2 or the metastables (zV^Re*. ( 2 D)N*, ( 3 P 2 Q )Ne*, ( 3 P 2 Q )Ar* 
(^D)O*, (a^4g)0 2 * can give performance (predicted by the theoretical 
specific impulse formula derived in Appendix A) superior to the maximum 
attainable for 0 2 -R 2 or an ^ ot ^ er conventional oxidizer-fuel combination. 

In particular, extraordinary theoretical performance is predicted for 
atomic hydrogen, the CH free radical, and the metastable excited states 
of helium, nitrogen, neon, and atomic oxygen with Hg. 

The essence of the manufacturing and storage problem is described in 
Section 5, prior to the detailed discussions of methods of production in 
Section 6 and matrix-isolation storage in Section 7. It is shorn in 
Appendix B that a very small characteristic time is associated with the 
heat pulse produced by a chance recombination reaction of matrix-isolated 
free radicals, and thus neighboring free radicals cannot absorb requisite 
energy from such a heat pulse for escape from their trapping sites. Be- 
cause the host matrix material will be vaporized along with product gases 
during the energy release, the matrix-isolation host species must also serve 
as the working fluid for a free radical or metastable propellant, and a 
decisive advantage is noted to accrue if the free radical or metastable can 
he produced in situ from one of the matrix host species H 2 » Ng, 0 2 , He, Ne 
or Ar. Deductive considerations show that among the propellant condidates 
with this property, atomic hydrogen is an H 2 matrix is most accessible to 
near-term development with existing technology. 

A review of current experimental work related to the manufacture of an 
atomic hydrogen propellant is presented in Section 8, where the two main 
technical aids now being employed to achieve higher concentrations of reactive 
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species with matrix-isolation, the use of very low temperatures and very 
strong magnetic fields, are discussed in detail. It is shown that im- 
portant progress toward the development of an atomic hydrogen propellant 
appears imminent with venerable methods of production and storage being 
augmented by these new technical aids. In Section 9 it is observed that 
the manufacture of an propellant, 15% atomic hydrogen by weight, may 

be possible by ultra-energetic hydrogen atom bombardment of solid 
(theory described in Appendix C) or by impregnation of solid E^ with radio- 
active phosphorus (theory sketched in Appendix D). The theoretical esti- 
mates presented in Section 10 show that matrix-isolated free radicals or 
metastables must be used directly without prior melting as solid propel- 
lants. Finally, the steady deflagration of an solid propellant, 15% 

atomic hydrogen by weight, is analyzed in Appendix E. 


2, Free Radicals 

Quantum mechanical analysis shows that a host of free radicals are 
likely to have adequate (ground-state) stability and sufficient specific 
energy of reaction to warrant further consideration as propellants 
(Schneiderman 170 , Kepford 107 ). The outstanding practical candidates 
for free radical propellants are shown in Table 1* in wmch the first eight 
entries , namely H, N f 0, CH, BH, MH, BHg and , 

are free radicals that have been studied experimentally for many years 
(Fimentel 157 , Milligan 134 ’ 137 , Windsor 210 , Ramsay 163 ), while the last 
three entries are the seemin^most promising additional free radicals that 
have been conjectured to exist by quantum mechanical analysis. 


’ J'TW 
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FREE 

RADICAL 

GROUND 

STATE 

CHEMICAL REACTION 

SPECIFIC ENERGY 
(kcal/gn) 

B 

2 S 

5 l/2 

i 

2H 

— h 2 

52.10 

N 

s 

S 3/2 


2N 

— N 2 

8.03 

0 

3 P 

i 

20 

— °2 

3.71 

CH 

V 

2CH 



2C (pftase) + H 2 

10.90 

BH 

/ 

/ 

3 n 

2BH 


2B + H 2 

6.19 

m 


2NR 


n 2 + h 2 

5.25 

bh 2 

h 

A i 

bh 2 


B + H, 

1 

5.18 

ch 2 

% 

ch 2 


c + h 2 

5.29 

He 2 0 

' V 

2He 2 0 

— 

4He + 0 2 

4.35 

HLi 

h* ~ 

2HL1 


H 2 + 2 LI 

5.46 

H 2 L1 

V 

U 

h 2 li 

— 

H 2 + LI 

3.94 


Table 1. Known and theoretically conjectured free radicals that 
may offer possibilities as propellants. 
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There is no experimental evidence that the latter entries , namely 
l^O, HLi and H 2 Li, can be produced in sufficient amounts to be 
useful as propellants (Kepford*^, Henneker^) , and since there is no 
advantage that would accrue in greater specific energy from the use of 
these latter entries, it is reasonable that attention in our survey should 
be concentrated on the three monatomic radicals H, N and 0, 
the three diatomic radicals CH, BH and 1TH, and 

the two triatomic radicals and CHg. 

3, Metastables 

An electronically excited state of an atom or molecule is said to 

be metastable if its radiative lifetime, is greater than a microsecond. 

Assuming that it is not possible to suppress the spontaneous electro- 

+ 

magnetic radiative decay of a metastable state , propulsion applications 

3 

would appear to require radiative lifetimes of the order 10 sec or 

i 

greater. A list of known metastables that offer even remote possibili- 
ties as propellants is given in Table 2. Known metastables with radiative 
lifetimes less than a millisecond or with relatively smaller values for 
their specific energy [ ■ (excitation energy of species) / (mass of species), 
expressed in column five of Table 2 by using the molecular weight values 

and the units conversion formulas 1 eV/partide a 23.0 kcal/ mole 3 

144 

have been excluded from a broader tabulation (Muschlltss ) . 


t It has been noted independently by J. Zmuid 2 inas and by the present writer that 
it may be possible to interfere destructively with the quantum mechanical 
probability amplitude for spontaneous electromagnetic radiative decay of an 
atomic metastable state, say by imposing a standing-wave electromagnetic field 
of critical amplitude and wavelength on metastables in a lattice. The theoretical 
analysis required to establish necessary critical conditions would be a fund- 
amental extension of the quantum mechanical calculations in Refs. 17, 50, 41-45, 
99, 111, 122, 214 and 215. 


t 


I 



i 


i 

l - 


l 


i 


4 

i 


i 




mental value exists 
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In the absence of a conceptual scheme for preserving the shorter-lived 

metastables, the electronically excited metastable species of practical 

3 2 3 

interest as possible propellants are 2 helium, D nitrogen, Pj ^ 

neon and argon, ^T> oxygen either with a nonreactive working fluid or 
with H 2 "after-burning”, and finally a 1 A oxygen with JU "after-burning" 

4c 

(i.e., the propellant 0 ^ — H^) . 

In addition to the basic metastfbles discussed above and shorn in 
Table 2, there are hydrides and other low molecular weight metastable 
compounds that can be made chemically by reacting (2 S-)He , ( J p )Ne , 

-L t j v 

( 0^ r *» ^ D)N*» or (^D)O* with ordinary (ground-state) hydrogen and 

other low molecular weight atoms* Such electronically excited metastable 
compounds have less specific energy than their metastable atomic progeni- 
tors, owing to the energy lost in the formation of the chemical bond and 
the greater molecular weight of the compound. Nevertheless, a metastable 
compound might be of considerable interest as a propellant if it were to 
have a radiative lifetime much greater than its metastable progenitor. 
Although one cannot prove a quantum mechanical theorem which demonstrates 
that this could never be the case, it is unlikely that compound-formation 

effects an increase in the radiative lifetime for spontaneous electro- 

+ 59 9 71 171 

magnetic decay of any metastable (Freed , Blxon , Guberman , Segal , 

t On the other hand, collisional reaction rates and decay rates 

associated with external electromagnetic perturbations can be 

suppressed by metastable compound-formation, as exemplified by the 

3 -5 3 + 

(2 S^)lle* states (with observed lifetimes about 1.5 x 10 sec) and (a T u ) 

He* states (with observed lifetimes greater then 0*1 sec) produced in 

electron-bombarded superfluid helium (see Refs. 40, 50, 75, 76, 81, 108, 

184, 185). The radiative lifetime of a metastable for spontaneous electro- 
magnetic decay provides a theoretical upper bound on the observed lifetime 
of a metastable that undergoes collisional reactions and/or Induced decay. 


j. 
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Roaenfield" 1 ^, Hcsbet^) . In fact, the interatomic Coulomb field felt 
by a metastable atom in a compound molecule usually increases the rate 
of electromagnetic radiative emission by a factor typically of the order 
10 , and hence the radiative lifetime of a metastable compound would 
ordinarily have a magnitude less than the radiative lifetime of its ato- 

_4 + 

mic metastable progenitor by a factor typically of the order 10 T . 

it 

Attempts to produce metastable helium hydride He H and, more recently, 

* * 

the metastable compounds OHe and ONe in significant amounts from 
(2 S^He and ( P 2 g)Ne for experimental applications have been un- 
successful (Kepford^^)tt Likewise, it has not been possible to produce 

and store electronically excited metastable states of free radicals in 

118 

significant amounts for detection (Linevsky ). Thus it appears most 
practical to concentrate future theoretical and experimental studies on 
the long-lived metastables of helium, nitrogen, oxygen, neon and argon 
that are listed in Table 2. 


♦ The interatomic Coulomb field felt by a metastable atom in a compound molecule 
generally has the order of magnitude of that due to a proton acting at a dis- 
tance of several Bohr radii t > oCm y*' Jn unit* fc»c* i * n <] eV4irs (is?) 4* 
This leads to the formal replacement of a fine-structure constant cc by a 
number of the order unity in the matrix-element for the transition, when one 
alters the vertex of a Feynman diagram for spontaneous radiative decay to give 
a Feynman diagram associated with the interatomic Coulomb field decay (see, 
e.g. » Ref. 50). Since the matrix-element is thereby increased by a factor of the 
order 0C" J , the radiative decay rate is increased by a factor of the order or* 
JO*" for the metastable atom in a compound molecule. 

^These metastables are theoretically auto-ionizing states (Ruseek 1 ^, 

Bhatia , Hiller , and Indeed experiments Indicate that auto-ionization 

Is concomitant with the formation of helium hydride (Neynaber 1 ^) , the 
dominant apparent process being He* + Hj, — +* HeH + H + (an energetic 

electron) with the resulting helium hydride ion unstable under subsequent 

electronic recombination. 
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4, Theoretical Perfoiuance 

For prescribed chamber conditions, nozzle geometry, end chemical 
composition of the exhaust gas, the theoretical performance of a pro- 
pellant is indicated by the associated value of the specific impulse. 

An approximate theoretical specific impulse formula for conditions typi- 
cal of those that would apply for an advanced propulsion system onboard 
a spacecraft (vacuum environment and essentially lossless isentropic 
equilibrium expansion through a nozzle with an expansion area ratio of 
60:1 and an effective constant specific heat ratio of 1.30) is derived 
in Appendix A with the result given by (A. 6), 

I ■ 265 (Q/l kcal/gm)* 5 sec, 
sp 

where Q denotes the overall net energy release per unit mass of propel- 
lant. Assuming that the free radical reactions in Table 1 go to comple- 
tion and the metastable species in Table 2 are quenched completely to their 
ground states, the energy release Q is given for any free radical or meta- 
stable propellant by multiplying the specific energy in Table 1 or Table 2 
by the weight fraction of the energized species in the propellant. Thus, 
for example, in the case of H with as the working fluid, an atomic hy- 
drogen weight fraction of .20 (meaning one free H atom for every two 

molecules) yields Q - 10.42 kcal/gm and I «* 857 sec, according to the 

®P 

approximate specific Impulse formula above. Table 3 displays Q and X 
values calculated in this manner for the best conventional chemical pro- 
pellants and the most promising free radical and metastable propellants. 

The weight fractions of the latter species are fixed at simple fractional values 
that give performance at least comparable to that of conventional propellants, and 
thus Table 3 shows the approximate weight fractions required for practical 
competitive use of free radical and metastable propellants. 


* 


FREE RADICAL 
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PROPELLANT 


WFIGET FRACTION OF 
OXIDIZER OR FREE 
RADICAL 


Q 

(kcal/grc) 


f o 2 - h 2 

Pj - h 2 
V»3-*2 


.825 

.938 

.889 


2.95 

3.16 

3.83 


/ H 


1.000 


52.10 


H 


( Inert working! 
fluid, e.g., J 
H, or He / 


.500 

.250 

.100 


26.05 

13.02 

5.21 


CH~( 

N-( 

, BH-( 

I NH-( 


) 

) 

) 

) 


b« 2 -( 

ch 2 -( 


) 

) 



^ 0 - H 2 


.500 

.500 

.500 

.500 

.500 

.500 

.500 

.500 


5.45 

4.02 

3.10 

2.62 

2.59 

2.64 

1.85 

3.64 


Table 3. Theoretical performance of the best conventional and the most 


I 

sp 

(sec) 


456 

471 

520 

1918 

1355 

959 

606 

620 

531 

467 

430 
427 

431 

361 

506 


promising 


free radical and metastable propellants 
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PROPELLANT 1 

WEIGHT FRACTION OP 

Q 

I 



METASTABLE 

(kcal/gm) 

sp 

(sec) 

+ » * / inert working \ 

l (2^S,)He -f fluid, e.g. , ) 

.500 

57.0 

2001 

.100 

11.4 

895 

1*2 

or Re / 

.050 

5.70 

633 

( 2 D)N*-( » 

) 

.500 

5.97 + 

648 


" ) 

.500 

.250 

9.52 

4.76 

820 

580 


" ) 

.500 

3.35 

486 

I - 4 . /inert, 

noncombus-^ 


3.26 TT 


} ( D)0 - I tible working 7 

.500 

480 

V fluid 

/ 




1 * 

CD)0 - H 2 


.500 

5.05 

597 

^ (a \ )0 2* " H 2 


.825 

3.53 

498 


Table 3 (continued). 


^Includes 4.02 kcal/gm resulting from the ground-state recombination 2>T — *» J? 2 
reaction that follows de-excitation. 

+ *Includes 1.85 kcal/gm resulting fro® the ground-state recombination 20 — *■ 0^ 
reaction but no contribution fro® a conventional oxidation reaction. 
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k 

Since I is proportional to and Q is directly proportional to the 
sp 

weight fraction of the free radical or metastab 1 species, it is a simple 
matter to figure the I value for any propellant in Table 3 with an al- 

sp 

tered weight fraction value. 

An inspection of Table 3 shows the extraordinary theoretical perfor- 
mance predicted for atomic hydrogen, the CH free radical, and the meta- 
stable excited states of helium, nitrogen, neon and atomic oxygen with Hj. 
Performance superior to the test conventional propellants is also indi- 
cated for atomic nitrogen with a weight fraction of .5 (meaning one N atom 
for each seven molecules if hydrogen is employed as the inert working 
fluid), while the free radicals BH, NR, and 0 would require weight 

fraction values in excess of .5 in order to surpass the performance of 
the best conventional propellants. Because of the relatively large mole- 
cular weights of (> 2 » ^ and 0^, an increase in specific impulse cannot 
be achieved by supplementing a free radical recombination that yields 

210 

as a product with a subsequent conventional combustion reaction (Windsor ) . 

On the other haud, a significant increase in specific impulse is attainable 

if the ground-state atomic oxygen recombination process, the T) metastable 

atomic oxygen quenching and recombination processes, or the a*A molecular 

6 

oxygen quenching process is supplemented by a conventional combustion re- 
action with Hj* as shown by the illustrative entries in Table 3. The para- 
mount issue is whether free radical and metastable propellants can be 
manufactured, stored, and used in practice in the required concentrations 
and amounts for rocket propulsion. 
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5. Manufacture and Storage 

Free radicals and metastables are produced by a large variety of 
physical and chemical molecular dissociation and electronic excitation 
reactions. Except under rather special conditions, however, such reactive 
species cannot be produced in high concentrations and/or in large amounts, 
as required for propellant manufacture. Moreover, for tl.e storage of 
highly reactive free radicals and metastables only one effective method 
is known, the so-called matrix-isolation technique in which the reactive 
species are trapped as isolated entities in ar> inert solid at a cryogenic 
temperature. Free radicals and metastables produced in a gas or liquid 
must be rapidly condensed and trapped by being frozen into normal or inter- 
stitial sites in an inert cryogenic lattice, while immediate storage of 
the species may be afforded by so-called in situ production methods which 
involve dissociation or excitation processes that take place exclusively 
within a pre-formed matrix. 

In the sections which follow, the methods for production of the re- 
active species and storage by the technique of matrix-isolation are re- 
viewed and considered from the standpoint of rocket propellant requirements. 
Despite the fact that no free radical or metastable propellant has been 
manufactured and stored in a practical quantity to date, the survey pre- 
sented here shows that important progress appears imminent with the venerable 
methods of production and storage being augmented by new technical aids. 
Indeed, the odds appear to favor the successful manufacture, storage, and 
use of an atomic hydrogen propellant within, the present decade. 

6. Methods of Production 

Table 4 provides a list of toethods that have been used to produce free 
radicals and metastables in gases, liquids, and solids. 


METHOD OF 
PRODUCTION 

REFERENCES FOR FREE 
RADICALS 

REFERENCES FOR 
METASTABLES 

PHOTOLYSIS: 

VISIBLE OR ULTRA- 
VIOLET RADIATION 

189, [157], [130— 138], 
[93 — 96], [164], [176], 
1, 47, 73, 79, 151, 158, 
159, 187 

149, [118], [125] 

X-RAYS 

[157] 

[125] 

tf-RAYS 

[157], [21], [202], 196 

1 

t 

ELECTROMAGNETIC 

DISCHARGES 

172, 11, 12, 16, 19, 51, 
56, 57, 146, 150, 26, 

80, 72, 200, 211 

145, 181, 118 

1 

I 

1 

THERMAL METHODS 

189, 25, 137, 179, 51, 
177, 18, 38, 205, 206 

121, 82, 3, 118 

CHARGE TRANSFER 
PROCESSES 

* 

145, 46, 201, 5, 
53, 178 

FLECTRON IMPACT 

[157], 203, 15 

[125], [118], 191-193, 

2, 126, 115, 29-31, 
183-185, 81, 83, 40, 

58, 13, 203, 15, 75, 

106, 108 145, 179, 209, 
204, 62, .3, 34, 69 

ATOM & MOLECULE 
BOMBARDMENT 

[157] 

l 

1 

86, 186, 82, [125] 


Table A. Methods of production of free radicals and metastables. Square brackets 
denote in situ production references, while all other references* describe 


production in a pas or liquid 
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Most useful for basic scientific investigations in which snail quantities 
of the reactive species are adequate, the very controllable method of 
photolysis utilizes visible or ultraviolet radiation produced by a high- 
intensity flash, an arc light, or a laser of suitable wavelength (belov 
the dissociation or close to the excitation wavelength of the species) . 

In photolysis, a so-called precursor species undergoes dissociation or 
excitation with the absorption of quanta in the intense radiation, there- 
by yielding the free radical or metastable (in the gas phase or in situ 
through an inert solid matrix at a cryogenic temperature, in usual appli- 
cations). The method of photolysis cannot be employed for the production 
of large amounts of free radicals or metastables, as required by propul- 
sion applications, because of severe limitations in the amount of energy 
that can be communicated to a medium by devices that emit visible or 
ultraviolet radiation. For example, a power level of the order of loV is 
required to energize reasonable amounts of solid to a useful P ro “ 

pellant mixture in a reasonable time (see Appendix C) ; such a required power 
level is several orders of magnitude above what can be delivered by a laser 
that emits radiation having a wavelength less than .276 ym (the threshold 
for dissociation of l^). 

Less developed than photolysis is the use of higher energy electro- 
magnetic quanta, x-rays and 2f-rays, which must effect dissociations and 
electronic excitations primarily through intermediate absorption processes. 
Although x-rays and Jf-rays allow for in situ production of the reactive 
species, the effective absorption cross-sections are small for light -atom 
solids and thus the higher energy electromagnetic quanta are not attractive 
energy sources for the manufacture of large amounts of free radical or 
metastable propellants. 
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Electromagnetic discharges (microwave, radiofrequency, etc.) and various 


thermal methods have been used widely for the gas phase production of free 
radicals and metastables. These methods are capable of producing free 
radicals and metastables at sufficiently high rates in the gas phase for 
possible use in the manufacture of practical amounts of propellants. In 
particular, free radicals can be produced copiously with glow discharges 
at low resultant temperatures which facilitate the required rapid conden- 
sation and trapping of the gas discharge in a cryogenic matrix. However, 
electromagnetic discharges and thermal methods are not applicable to 
in situ production and the more efficient concomitant trapping of the 

reactive species. Likewise, charge transfer processes, as exemplified 

• ^ ^ ^ 

by the helium ion react ion He + He — *■- (2 S^) He + He , have been found 
to produce metastables copiously under suitable conditions in the gas phase 
but not in a solid matrix. 

Most promising for the production of free radicals and metastables 
in the amounts required for propellants are electron impact techniques. 

The latter are effective for copious production in situ in a cryogenic 
matrix as well as in a cryogenic liquid, in a gas, or in (transverse bom- 
bardment of) atom and molecular beams. For in situ production of reactive 
species propellants by electron impact, it is necessary to generate free 
electrons in quantity throughout a matrix (see Sections 8,9 and Appendices 
.C, D), because even very energetic electrons have a very short range in a 
solid (e.g., an 'lectron with an inital kinetic energy of 1 MeV will pene- 
trate a distance of only about (O.A gm/cm )/p, where p is the density of 


the solid). 


Thlgh-e argy beam of electrons might be employed to energize thin wafers of a 
matrix -in a manner similar to the hydrogen atom bombardment in A PP e « d ^C,but 
power -delivery and charge-accumulation considerations favor neutral atom bom- 
bardment of a wafer. 


Methods for the production of reactive species in situ by atom and 

molecule bombardment appear to hold considerable promise in theory, but 

f 

etill have not been investigated with appropriate experimental research 
programs. At high energies for the atoms or molecules in the beam, se- 
condary electrons will be produced via ionization processes in the matrix 
and assist in the production of reactive species by electron impact re- 
actions (e.g., see Appendix C). 

All of the free radicals and metastables that appear in Table 3 have 

been produced in high concentrations by one or more of the methods in 

Table 4. Specific references for the efficient prjduction of H, N and 0 

27 

are given in the review article by Carstens et a3 . , for the production 

of CH, BH, NH, BH^ and CH 2 in the review articles by Milligan and 

Jacox^ 3 ^*^ 3 *^ 3 ^ and Ramsay*^, ar.d for the production of the metastables 

(2 3 S 1 )He* l ( 2 D)N*, ( 3 P 2 0 )Ne* and A <}v)0* and (a 1 A g )0 2 * in 

145 

the review article by Muschlltz . Moreover, the methods of production 
for the latter reactive species appear to be readily adaptable to scale-up 
and large-quantity manufacture. In fact, production efficiencies can be 
expected to improve with scale-up and the increase in volume-to-surface- 
area ratios, owing to the suppression of surface recombination and quenching 
reactions. Therefore, the essential problem is the stabilization and 
storage of the reactive species by matrix-isolation. 

~ ~ 

In 1960 Pimentel wrote (Ref. 157* p. 109): "The potentialities and limita- 

tions of atom bombardment are virtually unexplored .... The technique surely 
deserves continued exploration".. 
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7. Matrix-Isolation 

Important technical advances with the In situ and gas condensation 
applications of the matrix-isolation technique have been made in recent 
years by Pimentel , Woltner , Milligan , and Jacox 

Since the host matrix material will be vaporized along with reaction pro- 
duct gase6 during the energy release (see Section 10), the host species 
must also serve as the working fluid for a free radical or metastable 
propellant. Substances which have been used effectively for matrix-isola- 
tion of highly reactive species and have low molecular weights are Hj,^, 
Oj and the lighter inert gases He, !!e, Ar. The latter candidates for the 
dual role of matrix host and working fluid appear in Table 5 with the 
temperature values required to maintain the solids at prescribed vapor 
pressures (data extracted from Ref. 125, p. 192). With the exception of 
helium, all of these substances are solids that have molecular concen- 
trations close to .045 moles/cm^ at temperatures less than 4°K. 

Notwithstanding intensive development and applications of the techni- 
que, it has not been possible to secure a reactive species concentration 
greater than about 1% molar^ in a matrix by the now venerable in situ and 
gas condensation methods of trapping with temperatures of the solid as low 
as 4.2°K (Fontana*^* Wall^^, Brown^\ Windsor^**, Fite^, Bracknann^, 
Linevsky*^) . 


+ 

The molar concentration of about 1 % was obtained for the free radical NCN 

130-132 

in an argon matrix by Milligan with in situ flash photolysis of the 

precursor species cyanogen aside, K^CN. Also see Ref. 157, p. 106. 



TEMPERATURE (°K) REQUIRED TO MAINTAIN SOLID IN 
EQUILIBRIUM AT INDICATED "ALUE OF VAPOR PRESSURE 
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A principal difficulty that prevents the achievement of higher free 
radical and/or netastable concentrations is the imperfect isolation of 
free radicals in the direct deposition of a gas mixture or the imperfect 
isolation of precursor molecules for in situ production (Pimentel 161 ) + . 
Another reason why higher concentrations of reactive species have not 
been secured with direct deposition methods is that solids formed from a 
gas by rapid condensation at a very low temperature (i.e., less than 1/2 
the freezing point of the substance) have defect-laden crystalline or amor- 
phous structures which favor the loss of reactive species by diffusion 
away from trapping sites (Pimentel ) and inhibit the cooling of addi- 
tional condensed layers because of relatively poor thermal conductivity 
211 212 

(Windsor * ). A decisive advantage accrues, therefore, if the free 

radical or metastable can be produced in situ from one of the matrix 
host species in Table 5 after a nearly perfect lattice structure has been 
formed by gradual cooling to a very low temperature. Moreover, the 
exceptionally high thermal diffusivity of such an ultra-cold lattice makes 
it possible to control the thermal side-effects of in situ production. 


A systematic experimental investigation aimed at identifying precursors 
which Isolate themselves preferentiall y (prior to in situ reactive species 
production) has not been performed to date, however. 

'The critical temperatures for the escape of atoms from traps in a rapidly 
formed (imperfect) lattice generally lie in a band from about 1/7 to 1/2 of 
the substance* 8 freezing point, thus from about 9°K to 32°K for N in 
Ng and from about 2°K to 7°K for H in Hj (see kef. 157, pp, 75-78). 

For a nearly perfect lattice formed by gradual cooling of the substance, some- 
what higher values for the bands of critical temperatures are indicated. 
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which might otherwise nullify the trapping effectiveness of a matrix with 

local heating that promotes the escape of reactive species from trapping 

sites. Those free radical and netastable propellant candidates in Table 

3 which can be produced in situ from matrix host species in Table 5 are 

2 

atomic ground-state hydrogen, atomic ground-state and D metastable 

nitrogen, atomic ground-state and *D metastable oxygen, metastable 

3 1 

helium, P, n metastable neon and argon, and molecular a A oxygen. Since 

practical competitive weight fractions (see Table 3) for ground-state 

2 

or D metastable atomic nitrogen are unlikely to admit stabilization in 

an N 2 matrix under the most favorable conditions (see Section 8) , ground- 
2 

state and D metastable atomic nitrogen might be precluded as a prime 
candidate. The requirement of a radiative lifetime greater than 2 hours 
would preclude the metastables of atomic and molecular oxygen and very 
probably the metastables of neon and argon as well^. Finally, the severe 
technical difficulties involved in maintaining helium aa a solid (see 
comment in caption under Table 5), and perhaps also the limited availa- 
bility and high cost of helium ( $35/kgm) might preclude the helium 

metastable. Our prime candidate for a propellant is assuredly ground-stete 
atomic hydrogen in an H 2 matrix. 

Particularly noteworthy, the very low concentrations of reactive 

species obtained experimentally with venerable matrix-isolation techniques 
are well below what is predicted to be attainable by theoretical analysis 
that assumes an ideal lattice and the random static isolation of free radi- 
cals, subject to the proviso that two adjacent (nearest-neighbor) trapping 

" " " 11 1 " 1 ' ~ 

The observed lifetime of a matrix-isolated metastable species can be con-, 
slderably leg® than the (vacuum-isolated) radiative lifetime because of 
the perturbing Influence of the lattice (see Ref, 89, pp. 329-334). 


! 
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sites cannot both accommodate reactive species (Jackson 89 * 90 , Golden 70 , 

32 

Chessin ). According to such theoretical investigations of the problem, 

mole fractions in the neighborhood of 15% should be attainable with a 

nearly perfect lattice structure, in situ production, and ideal isolation 

(that admits only nearest-neighbor recombinations). However, 

91 

Jackson has given a dynamical stability analysis based on a phenomeno- 
logical continuum model which appears to indicate that local 
heat diffusion from a chance reaction of a pair of free radicals (e.g., 

H + H — 2 or N + N — ^NjJmay limit stable matrix-isolated concentrations 
of reactive species to less than about 2 or 3% molar. It is shown in 

Appendix B that Jackson's model leads to a characteristic heating time of 

—14 Q 

the order 10 sec if the matrix is maintained at 4 K or at a lower 

(realizable) temperature. Clearly, free radicals in the vicinity of a 
reacting pair cannot absorb requisite energy and escape their trapping 
sites during the very small characteristic time associated with a heat pulse 
in a matrix at 4°K or at a lower temperature. A quantum mechanical treat- 
ment of the trapping stability problem, which takes account of the discrete 
character of particle and energy transport processes, would be of considera- 
ble Interest and importance. 


8. Current Experimental Work 


To achieve the higher concentrations of free radicals and/or metaatablea 
required for useful propellants, it appears necessary to augment the vene- 
rable matrix-isolation method with sunplementary technical aids that assist 
In the stabilization of the reactive species. Two main supplementary tech- 
nical aids are now being employed in experiments, namely, the use of very . 
low temperatures (in the range 0.1°K to 1.5°K) produced by advanced cryogenic 
apparatus and the use of very strong magnetic fields (in the range 50k<5 to 
lOOkO) produced by superconducting magnets. At temperatures below 1.5°K, 
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one can expect a dramatic improvement in the trapping effectiveness of 
a matrix compered to what it is at 4.2°K, the liquid helium temperature 
of earlier matrix-isolation experiments. In particular, for hydrogen 
atoms in a nearly perfect lattice, the lower edge of the band of cri- 
tical temperatures for the escape of atoms from traps is somewhat greater 
than 2°K, which implies that the trapping processes should proceed 
favorably at temperatures below 1.5°K. Very strong magnetic fields can 
saturate the alignment and make parallel the spins of unpaired electrons 
in free radicals, thereby inhibiting radical-radical recombination re- 
actions. 

The use of strong magnetic fields to assist in the stabilization of 

matrix-isolated free radicals is not a new idea (see, for example, 

_ 104 . 210,211 51 „ 16 v , . , t 

Jones , Windsor , Fite , Brackmann ) , but only recently has the 

cryogenic and superconducting magnet technology been available for the 
attainment of conditions favorable for the stabilization of atomic hydro- 
gen, according to quantum mechanical theory. Table 6 shows the low tempera- 
ture-strong magnetic field combinations predicted theoretically (Ref. 104, 
p. 25) to yield substantial saturation of electron spin-alignment in ato- 
mic hydrogen and thus to inhibit the recombination reaction 2H — *- 
(with only a total electron spin zero bound-state admissible for molecular 
hydrogen). However, the combinations of low temperatures and strong mag- 
netic fields In Table 6 are neither necessary nor sufficient, for atomic 

For free hydrogen atoms at temperature T in a strong magnetic field of 
magnitude B, the recombination rate Is suppressed by the Boltzmann pro- 
bability factor for non-alignment of the electron spin in a hydrogen atom, 
given by e ” 2 vB/kT ^ « 1°K/15 kO Is the spin-magnetic moment 

of an electron divided by Boltzmann’s constant* The values in Table 6 
obtain if the probability factor is required to equal .25 x 10**, a number 

that makes the lifetime of the atomic hydrogen sufficiently long for 
practical utilization (Ref. 104). 



Table 6. Combinations of low temperatures and strong magnetic fields pre- 
dicted theoretically to stabilize atomic hydrogen by nearly con- 
plete alignment of electron spins. 

hydrogen trapping in an matrix will assist stabilization, while the 

small interaction between the magnetic dipoles associated with the elec- 

104 

tron and proton spins in a hydrogen atom (Jones ) and electron-spin 

139 

multipole radiative transitions (Mizushima ‘ ) may promote recombinations 
by electron spin-flips in the course of interactions between atoms with 
electron 6pins parallel initially. The rate of the latter phenomenon 
is quite difficult to estimate for electron spin-alj-gned hydrogen atoms in 
an Hjj matrix* . 

Exciting positive experimental results for atomic hydrogen stabiliza- 
tion have been reported recently. In a painstaking research program, 

80 

Hess has produced significant concentrations (apparently about 10£ molar) 
of atomic hydrogen in the malrlx-isolated condensate from a hydrogen glow 
discharge cooled very rapidly to about 1.3°K in a magnetic field of appro- 
ximately 50kG . 

The required quantum mechanical calculation would have features in common 
with the analysis reported in Ref. 156. 

**Eess®® has conjectured that large clusters of spin-aligned atomic hydrogen 
may have formed and precipitated out in the condensate as atomic hydrogen 
domains, but this seems unlikely in view of recent theoretical estimate made 
by Etters^®. 


i 
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A copy of the Eess apparatus has been constructed at NASA's Lewis Pesearch 
Center and experiments now in progress will clarify the effectiveness of 
this method for atomic hydrogen propellant manufacture and storage. In a 
thermal sensing experiment, a mixture of H and Hj (or dissociated and 
molecular deuterium) is deposited on a nearly thermally isolated collector 
at about 1.2°K in a magnetic field of ICOkG; the recombination of atomic 
hydrogen is induced by electrical heating and the energy release is mea- 
sured by a carbon thermometer. In a mass spectrometer sensing experiment, 
the bee.'-, of particles from either the H-H^ source or the collection re- 
gion is chopped by a toothed wheel and phase-sensitive detection is possi- 
ble with an in-line-of-sight quadrupole mass spectrometer which is capable 

of recording particles with molecular weight 1, Pesults from t’ sse experi- 

22 

ments are expected to be available later this year (Brown ) . 

Work in progress will also determine the effectives js of internal in 
situ free radical production at ultra-low temperatures in a strong magnetic 
field. Here the idea is to use the $-ray electrons emitted in the decay of 
tritium to produce atomic hydrogen by electron impact. In experiments now 
being performed at Lewis Research Center, ordinary (molecular weight 2) 
and tritium molecules (weight 6, half-life 12.26 years, mean energy of 
8-ray electrons 5.7 keV) are frozen together in a dilution refrigerator 
at Q.1°K and a 35kG magnetic field is applied to the sample, thus providing 
very favorable conditions for stabilizing H atoms that result principally 
from secondary electron impact dissociations of molecules in the matrix 
(see Appendix D). Measurement of the atomic hydrogen concentration is 
accomplished by determining the magnetic susceptibility of the solid. This 
experiment will also yield results presently (Brown ). 
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9. Alternative Propos als fo r Ma n ufactur in£_ H-H a Prop ellant s 

In addition to the methods of manufacture and storage of an atomic 

hydrogen propellant that are under current investigation, there are tvc 

alternative schemes vhich also pppear to warrant experimental study. The 

first involves the use of a cyclotron to produce a 10 MeV hydrogen atom 

beam which bombards solid at a temperature below 1°K in a very strong 

magnetic field. This possible method for the manufacture of an F-II P ro ~ 

pellant (about 15% atomic hydrogen by weight) has been worked out in a 

preliminary theoretical fashion by the present writer and is described 

here in Appendix C. The second scheme involves the use of the radio- 
32 

nuclide P (with more than 300 times the specific activity of tritium) 
to produce an K-llj propellant through secondary electron impact dissocia- 
tion reactions. Secondary electrons will be proliferated by 6-ray elec- 
trons (mean energy 690 keV) in an matrix impregnated v:ith .002* molar 
32 

P atoms, and conditions trill be favorable for the trapping of free hydro- 
gen atoms that result from dissociations if the solid Is maintained below 
1°K and a strong magnetic field is applied during the radioactive "curing" 
of the propellant. Essentially an extension of the tritium IT— propellant 
concept to higher energetics and production rates, this possible method 
for the manufacture of a H-Hj propellant (at least 15% atomic hydrogen by 
weight) has been sketched by the present writer in Appendix D. Both of 
these two schemes offer promise of providing an atomic hydrogen free radical 
propellant that v?ill yield a specific Impulse of about 740 sec for rocket 
engine conditions that make formula (A. 6) applicable. 
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1 0. Energy Relea.s • 

For all of the highly energetic free radical recor'.bination reactions 

in Table 1 the activation energies are close to zero, and these reconbi- 

nation reactions proceed rapidly even at cryogenic temperatures if the 

174 1 

species are mobile in the gas or liquid phase (Sinha , Adrian ) . Accord- 
ing to elementary collision reaction theory, the rate of such a free radi- 
cal recombination reaction is given by f a er n v where c is the velocity- 
averaged total cross-section for the reaction, n is the concentration of 

the free radicals, and v is their mean velocity in the gas or liquid. Now 

— -18 2. 
the cross-section a is certainly greater than 10 cm for the radical- 

radical reactions of interest at low energies^ and it is required that n 
be greater than 10 free radicalsXm' for a useful propellant in the stored 
(matrix-isolated) cryogenic solid state. Thus, if a cryogenic solid con- 
taining trapped free radicals were to melt and attain a temperature of, say, 

o — 4 

30 K and an associated mean particle velocity v > 10 cm/sec, then the 

characteristic time for the free radical recombination reaction would be 

f ^»s= (ff n v) ^ < 10 ^ sec. Quite clearly, in a time less than a tenth of 

a microsecond it would not be possible to pump and utilize the free radical 

mixture as a liquid propellant. Hence, it is necessary to use matrix- 

isolated free radicals directly without prior melting as solid propellants. 

The necessity of using matrix-isolated metastables as solid propellants 

can be established by an argument similar to the one given above for free 

radicals. Electronically excited metastable species can react chemically or 

7 

Two orders of magnitude larger than the conservative general lower bound 

evoked here, the total cross-section for the reaction H + H — H_ is about 

—16 2 ^ 

1.1 x 10 cm at thermal energies (< 0.5 eV), with account taken of the 

fact that the reaction proceeds only in the total electron spin zero channel. 


i ip < 
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undergo an inelastic collision with virtually any atom or molecule, in- 
cluding the inert gases in their ground states. Three types of colli- 
slonal quenching or de-excitation processes are possible: (1) collisions! 

processes leading to thermal (increased molecular translational energy) 
or radiative release of the metastable species excitation energy (see Refs. 
23, 24, 49, 68, 105, 110, 114, 124, 142, 144, 145) , (2) ionization reactions (see Refs. 
4,20,35,85, 87, $8,113,128,144,154,162,180,188), and (3) excitation 

transfer reactions (see Refs. 68,112,166,173). The cross-sections for all 

-15 2 

of these reactions are large, typically of the order 10 era at low to 
moderate particle velocities. Although the latter ionization and excita- 
tion transfer reactions do not in themselves effect the release of exci- 
tation energy in thermal or radiative form, the products of the latter 
reactions quickly undergo subsequent electronic recombination and energy- 
liberating collisional processes at the high species concentrations of 
Interest. Thus, if a cryogenic solid containing a useful concentration 
of trapped metastables were to melt and attain a temperature of, say, 30°K, 
the characteristic time for overall effective quenching and sensible 
energy release would be of the order or less than about a nanosecond, pre- 
cluding use of a cryogenic solid containing metastables except directly as 
a solid propellant. 

Ignition of a free radical or metastable solid propellant can be accom- 
plished by any means (electrical or other) that warms a surface layer of 
the cryogenic matrix to a temperature at which the rate of heat release 
from recombination or quenching of the reactive species exceeds the rate of 
heat 1068 by thermal diffusion. Generally, this will be the case if a 
surface layer of millimeter thickness is mainlined above the critical tem- 
perature for significant escape from trapping sites (expected to be in the 
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vicinity of 4 K for hydrogen atoms in a well-formed H matrix) for a 
tine of the order of a millisecond. Because retastable species can also 
be quenched efficiently by suitable static electric fields or by elec- 
tromagnetic radiation of appropriate wavelength (see Refs. 17,84,99,111,117, 
141,155,214,215), more controllable nonthermal ignition systems are con- 
ceivable as alternatives for metastable solid propellants. 

Under proper conditions, the deflagration that follows ipnition of a 

free radical or metas table propellant can be expected to be steady and 

stable, albeit very rapid, with the recombination or quenching reaction 

taking place almost exclusively in the fluid (very light liquid and gas) 

adjacent to the melting surface of the solid matrix. An elenental volume 

of the solid will be warmed to the melting point in a tine much less than 

a microsecond, a duration too brief for reactive species to absorb energy, 

escape from their trapping sites, diffuse through the solid and undergo 

reaction. Heat transport into the interiors of the matrix will be restrained 

because the solids of interest have thermal diffusivities which decrease 

o 

rapidly with increasing temperatures above 1 K. For example, the thermal 

diffusivity of solid hydrogen is given approximately by (Ref. 36, p.13) 

D & 1.5 x 10^ x (exp - T/.845°K)cm^/eec in the temperature range 4°K < T J 
o 2 2 

12 K and remains of the order l(r cm /sec (with roughly proportional decrea- 
sing thermal conductivity and specific heat) as the temperature decreases 
from about 4°K to below 1°K, Moreover, because the solids of interest are 
generally transparent to the bulk of electromagnetic radiation in the Infrared 
and visible wavelengths, the absorption of radiation emitted by hot gas in 
the chamber Is not expected to produce a significant temperature rise in the 
depths of the solids. 


4 
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Straight- forv&rd order-of-mngnitude considerations show that the 
burning velocity of a free radical or netastable solid propellant (i.e., 
the rate of linear regression of the melting surface) will be determined 
by the heat flux from the hot fluid to the surface via the equation 



where \ Is the thermal conductivity of the liquid at the melting point 

(surface temperature) of the solid, AT is the temperature rise in the 

fluid in a small distance Ax (of the order of the local mean free path 

for molecules in the fluid) measured from the surface, V is the burning 

8 

velocity of the solid, p is the density of the solid, and h, is the en~ 
thalpy required to warn and melt a unit mass of the cold solid. The defla- 
gration of an solid propellant (15% atomic hydrogen hy weight) is dis- 

cussed quantitatively in Appendix E, where the steady-state burning velocity 
V * 2.1 m/sec is predicted for the propellant at a chamber pressure of 

9 

p - 100 psla. As shown by the analysis in Appendix E, the heat flux to the 
solid is moderated by the large velocity of the flow and the rapid expan- 
sion of the fluid with heat release from the recombination reaction. 

In the handling of a free radical or meta3tahle propellant, careful 
precautions must be taken to Insure that a sufficiently lovr stabilizing tem- 
perature (augmented by a sufficiently strong magnetic field, in the case of 
a free radical propellant) is maintained uniformly throughout the matrix at 
all times. It is also clear that contact of the solid propellant with any 
material which might catalyze energy release or absorb the reactive species 
from the matrix must be avoided. For example, It Is known that platinum, 
nickel and certain other metals can catalyze the recombination of atomic 
hydrogen at low temperatures. Moreover, atomic hydrogen is absorbed readily 
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by and diffuses interstltially through virtually all ferrous alloys In 

175 88 

contact with molecular hydrogen (Smialovski , Interrante Ivith the dis- 
sociation H 2 — *-2H catalyzed at any hydrogen-steel interface, and it may 
not be possible to contain an ultra-cold H-H^ solid propellant in a high- 
strength steel casing. Solutions to the handling problems associated with 
free radical and metastable solid propellants should be obtainable by per- 
forming systematic experimental research with sample quantities of the 
cryogenic solid propellants. 


11. Recommendations for Future Work 


The following problems are of prime importance and worthy of detailed 
theoretical analysis in the near future: 

CL> A quantum mechanical treatment of the trapping stability problem 

which takes account of the discrete character of particle and energy 
transport processes in a matrix. 

(2) The specialization of (1) for atomic hydrogen in an matrix, extended 
to include the quantum mechanical effects of a superimposed strong 
magnetic field on the stabilization at very low temperatures. 

(3) The extension of (2) to include the effects of heat transport into a 
pre-stabilized H-H 2 propellant that undergoes surface deflagration. 

It is recommended that the current atomic hydrogen experimental program at 
LRC be expanded to include other basic experiments which will provide data for 
testing the theoretical results of (2) and (3). 
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Propellants composed of any of the free radicals H, II, 0, CH, PH, 

KH, BH 2 , CH 2 or the metastables (2 3 Sj)He*, ( 2 D)>I*, ( 3 P 2 0 )I T e*, ( 3 ? 2 Q )Ar*. 

(^D)O*, (a^A^)0 2 * can, In theory, give specific impulses well above the 
maximum attainable for 0 2 - H 2 or any other conventional oxidizer-fuel com- 
bination. In particular, extraordinary theoretical performance is pre- 
dicted for atomic hydrogen, the CH free radical, and the metastahle excited 
states of helium, nitrogen, neon, and atomic oxygen with H 2 * 

A large variety of molecular dissociation and electronic excitation 
reactions can be employed to produce free radicals and tr.etastables in gases, 
liquids and solids, but the reactive species cannot be produced in high 
concentrations and/or in large amounts except under rather special conditions. 
Moreover, there is only one method known to be effective for the storage of 
these reactive species, nanely, matrix-isolation. In the latter technique, 
the reactive species are trapped as isolated entities in an inert solid at 
a cryogenic temperature, either by rapid condensation of a gas mixture or 
with in situ production in a pre-formed matrix by electromagnetic quanta or 
energetic particle bombardment. Although no free radical or metastable pro- 
pellant has been manufactured and stored in a practical quantity to date, 
Important progress appears imminent with the venerable methods of production 
end storage being augmented by new technical aids. 

Electromagnetic discharges and various thermal methods are 
capable of producing free radicals and metastables at sufficiently high 
rates in the gas phase for possible use in the manufacture of practical 
amounts of propellants, but these methods are not applicable to in situ 
production and the more efficient concomitant trapping of the reactive 
species. Particularly effective for the production of free 
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radicals and netastables in situ , electron inpact techniques are rest pro- 
mising for the manufacture of reactive species propellants. The short 
range of energetic electrons in a solid requires that they he generated in 
quantity throughout the matrix. Methods for the production of storable 
reactive species In situ by atom and molecule bombardment also hold consi- 
derable promise and warrant future experimental study. 

Since the host matrix material will be vaporized along with product 
gasea during the energy release, the host species must also serve as the 
working fluid for a free radical or metastable propellant. Candidates for 
the dual role of matrix host and working fluid are ^ 2 * ^ e ' *' e ari< ^ 

Ar. A decisive advantage accrues if the free radical or raetastable can be 
produced In situ from one of the latter matrix host species. Among the. 
propellant candidates with this property, atomic hydrogen in an matrix 
Is most accessible to near-term development with existing technology. 

Two main supplementary tech ileal aids are now being employed in experi- 
ments to achieve higher concentrations of reactive species with matrix- 
isolation, namely, the use of very low temperatures (in the range 0.1°K to 
1.5°K) produced by advanced cryogenic apparatus and the use of very strong 
magnetic fields (in the range 50 kG to 100 kG) produced by superconducting 
magnets. A significant Improvement in the trapping effectiveness of a 
matrix is expected at temperatures below 1.5°K, while very strong magnetic 
fields can prevent radical-radical recombination reactions by effecting 
the parallel alignment of unpaired electron spins in free radicals. Signi- 
ficant concentrations of 'atomic hydrogen, apparently about 10% molar, have 
been obtained in the matrix-isolated condensate from a hydrogen glow dis- 
charge cooled very rapidly to about 1.3°K in a magnetic field of approxi- 
mately 50 kG. Similar experiments in progress with magnetic fields as 
largo as 100 kG vill clarify the effectiveness of this method for atc"ic 
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hydrogen propellant manufacture and storage. The use of a radioactive 

32 

p~ray emitter isotope (tritiun under current experimental study, P as 

a proposed alternative) may yield F*-H 9 propellants (with I =• 740 sec 

32 

estimated In the case of P ) by secondary electron inpact dissociations 
of P *2 In an Impregnated matrix maintained below 1°K (prospectively at 
0.1°K in the case of tritium) in a strong magnetic field. Another method 
for manufacturing an H-Hj propellant ( — ' 15% atomic hydrogen by weight , 

I & 740 6ec) involves bombardment of with a cyclotron-produced beam 

8p L 

of 10 MeV hydrogen atoms. 

Estimates of the characteristic times for generic free radical recom- 
bination reactions and generic metastable quenching reactions show that 
the matrix-isolated reactive species must be used directly without prior 
melting as solid propellants. Under suitable conditions, the deflagra- 
tion of a free radical or metastable propellant is expected to be steady 
and stable, albeit very rapid. A burning velocity of about 2.1 m/sec is 
predicted for an H-Hj propellant (•*> 15% atomic hydrogen by weight) at a 
chamber pressure of 100 psla. The handling problems associated with such 
solid propellants are believed to be amenable to experimental solutions. 
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Appendix A: Theor et i c *l Specific Impulse kor nula 

Let Q denote the overall net energy release per unit mass of propel- 
lant end n denote the overall efficiency of the rocket engine in conver- 
ting this energy release to kinetic energy of the gas at the exhaust plane 
of the nozzle. Assuming that the sensible enthalpy of the cold propel- 
lants is negligible^ compared to Q, the specific Inpulse is given by 

I - / g, (A.i) 

op 

2 

where g = 980 cm/sec is the acceleration of gravity at sea level. The 
overall efficiency for energy conversion in (A.I) can be expressed as 

n * n j . i n, (a. 2) 

ideal losses 

in which ^ • 1 — (temperature of gas at the exhaust plane of the 

nozzle) / (chanber temperature of the gas) is the idea] thermodynamic effi- 
ciency of the rocket engine, and the factor nj 08£ , e8 ($l) accounts for real 
fluid dissipative effects. Ve have 

■W + (A.5) 

according to the theory for inviscid, iscntropic, one-dimensional flow of 
a perfect gas with # denoting the effective constant ratio of specific 
heats and M denoting the ideal Mach number of the gas at the exhaust 

C* - 

plane of the nozzle. With optimum design the practical efficiency factor 
n losses * n ^*2)» associated with possible reductions in the exhaust velo- 
city and * f ach number due to unrecoverable hert losses to the chamber and 

o 

nozzle walls, real gas viscous flow effects, etc., can he brought to 


tFor example, the sensible enthalpies of H«> and 0^ at their normal boiling 

points (20°r. and *-0 o li, respectively) are both less than .03 keal/gm and 
hence negligible compared to their atoichionetrlc energy release Q ■ 3.18 
keal/gm. 


I i 
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a value close to unity . From inviscid, lsentropic, one-din ensional 
flow theory we also obtain the relation between the e:<pansion arm ratio 
of the nozzle, denoted here by R, and the ideal Mach number of the pas 
at the exhaust plane, 

R - M [2(b+l) _1 + (*-l) (»41) -1 M * (A. 4) 

cx ex 

which can be inverted for R»3 to yield 

c*. s) 

Thus, for a spacecraft with an expansion area ratio R ** 60 and e. .haus t 

gas with an effective constant ratio of specific heats X «= l.?5, we 

find that (A. 5) and (A. 3) yield M a ■ 2<t.O and n., , * .750, while for 

ex ideal 

R ■ 60 and 1.30 (A. 5) and (A. 3) yield * 27.9 and n i(j ea i * .808 

respectively. If we assume the latter values for R and X and set the 
efficiency factor 0j ogses ®<jual to unity in (A. 2), we obtain the theoreti- 
cal specific impulse formula from (A.l) and (A. 2), 


^Provided that thermodynamic equilibrium is maintained between the molecular 
(internal and translational) degrees of freedom during the expansion of the 
gas through the nozzle, energy partitioned Into the vibrational, rotational 
and electronic excitations la taken Into account In the value for X . Since 
X decreases tovard unity as the number of affective Internal degrees of 
freedom lncreasea, the Ideal thermodynamic efficiency la reduced by en 
Increase In the number of Internal degreea of freedom according to (A. 3) and 
(A. 5) with R fixed. On the other hand, energy losses due to nonequlllbrlun 
finite- time relaxation of the Internal degreaa of freedom must be taken 
into account through the factor nj 08geg . 
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I = 265 (Q/l kcal/gm)* 3 sec, 

Bp F 

[for R - 60. *. 1.30, o lo5seB - l]. 


(A. 6) 


Formula (A. 6) is used to estimate the performance of propellants in the 
body of this report. For example. In the case of the conventional chemi- 
cal propellant hydrogen-oxygen system with the optimum mixture ratio 4.7 
of 0, to H* by weight, we have Q - 2.95 kcal/gm and (A. 6) gives I » 456 

* t 8p 

sec. 


Appendix B: Local Heat Pulse Produced by a Chance 

Reaction of Free Radicals in a Matrix 

91 

Following Jackson , let W denote the energy transmitted tc the lattice 

by the recombination reaction of a pair of free radicals. Since only a 

fraction of the total energy release from such a recombination reaction will 

be absorbed by the lattice, we have W < 4.5 eV for R + and W < 9.8 

eV for N + N-fr-Ng’ bet denote the minimum local steady-state temperature 

required for freeing a radical from a trapping site and T denote the Ini- 

o 

tlal uniform temperature of the solid matrix. Then trapped radicals may 
possibly be freed by the approximately Gaussian heat pulse that emanates 
from a chance reaction through a spherical region of radius 

r £ sk [w/nc <T f - T o )] 1/3 (B.l) 

3 

In which nm.045 moles/cm i_. the concentration of molecules (e.g., Hj or 
H 2 ) In the solid and c a [(T f + T o )/20°x] 3 cal/m©le - °K • 2.6 x 10 19 C(T f + T o ) 

/20°X] 3 eV/mole - °K is a simple approximation for the mean molar specific heat 
of the solid (proportional to T 3 at temperatures below about 15 X) In terms of 
the mean temperature £(lf + T C ) over the thermal range of Interest (see Ref. 91). 
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(B • 2) 


(B.3) 
(B.4) 

It Is evident that radicals in the vicinity of a chance recombination 
reaction cannot absorb requisite energy and escape their trapping sites 
during the small characteristic time of order 10”*^ sec in (B.3' and (B.4), 

for a radical would have to acquire a velocity of the order 10** cm/sec to move 

an Angstrbm during the characteristic time of the heat pulse. Hence, a quantum 

mechanical treatment of the trapping stability problem is required. 

Appendix C: Manufacture of an H-Ha Solid Propellant 

by Pltra-Fnergetic Hydrogen Atom Bonbardment 

A conventional cyclotron of 1950 vintage, with e magnet of 100 inch 

diameter end redlofrequency power supplied at about 250 kV, will deliver 

an external (extracted) beam of lOHeV protons at a steady current of about 

•2 mA * 1.2 x 10** protons/scc (implying a beam power of about 192 kV) . 

Such a beam of high energy protone can be made to have a uniform flux over 

3 2 

a cross-sectional area of about 3 x 1C cm (by eleetrle or magnetic field 

■ 

dispersal focusing) and to pass from vactxaa through a chamber filled with 


Letting D denote the thermal dlffusivity of the solid, it follows that 
trapped radicals within ii stance r^ of the recombination reaction experi- 
ence the heat pulse during a time duration given approximately by 

2 /v\ ^ T»i / f*r> - o?/3, 


t f s r f VD & [v/nc(T f - T o >] /d. 


For the typical parameter values W * 2 eV, ■ 10°K, D ==• 150 cm^/sec at 

T & 4°K and D & 100 cm^/sec at T £ 1°*’, ve obtain 
o o 


r^ & .95 x 10 ^cm ^ 

t^ & .60 x 10 *^eec ** 
-6 

a 1.05 x 10 cm T 
t* at 1.10 x 10~*^sec ' 


at 


at 


T * 4 K 
o 


T & 1°K. 
o 
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17 3 

cesiun vapor at a density of al out 10 cosiuii cic: -/cr . Since the 
total cross-section for the charge-exchange process leading to anv neu- 
tral (possibly electronically excited) state. o r hydrogen, H + + fe — •- 
B + Ce* t can be estimated for 10 'AeV protons by exti a; c let ion of lower 
energy data (Spiess ) to be at least 10 ’ cm , It fo] lot’s that almost 

*11 of the protons will be converted to neutral hydrogen atoms If the 

cesium vapor neutralizing chamber has a length along the bear of about 

86 

10 cm (see Hughes for a detailed discussion of the technique). In this 

manner it is possible to produce a beam of 10 MeV hydrogen atoms wi* - ’ - a 

13 7 

uniform flux of about 4 x 10 particles/cn -sec over a cross-sectional 

3 2 

area of about 3 x 10 cm . 


Now suppose that a 10 liter volume (equivalently, .73 kpn) of 1 quid 
1»2 is cooled gradually through its freezinp point and down to a temperature 
below 0.5°K to yield a wafer of so3id Fj about 3 x 3 O^cn^ in area and 2.7 cm 

in thickness. The wafer of solid H^, a rigid lattice composed of about 

26 

2.14 x 10 molecules, will have a nearly perfect crystalline structure, 

good thermal conductivity, very low specific heat, and thus very high thermal 

36 127 

diffusivity (Contreras , Miller ). 

Let us consider bombardment of the ultra-cold wafer of solid H 2 by the 
beam of 10 MeV hydrogen atoms described above. Having a velocity of about 
4.4 x 10^ n/sec, the 10 MeV hydrogen atoms will enter the aolld without the 
significant backseat ter and reflection observed experimentally at .nuch lower 
bombardment energies (Brackrann** , Fite**). This la because at the high 
velocities of atoms in the beam the cross-sections for F-Hj scattering and 


reaction procet 


»ry small (l.e., well below the geometrical 4 x 10* 


cm; and scattering amplitudes are sharply peaked in the forward direction. 
Because the Hj molecules in the solid lattice era separated by a distance 
(about 3*) which la roughly three times the diameter of an Hj molecule, a 
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■f* 

10 MeV hydrogen atom will travel about 2.0 cn in the solid before beinp 

brought to rest by energy-absorbing processes. Most of the energy of the 

10 MeV atoms will be absorbed by ionization reactions which take on the 

average 38 eV per ion-electron pair and yield energetic free electrons 

that move preferentially in the direction of the hydrogen atom beam. There 

will also be knock-on dissociation reactions II + II^ - *" 311 by the primary 

i-ydrogen atoms, yielding secondary energetic chain-branching hydrogen 

atoms which also move preferentially in the direction of the original beam. 

The secondary energetic free electrons and hydrogen atoms will also give 

up tbeir energy primarily through ionization reactions and dissociation 

+t 

reactions . In the course of these reactions, the excitation of lattice 
vibrations (associated with the sensible heat content and temperature of 
the solid) will not be favored quantum mechanically if the temperature is 
maintained below 1°K during bombardment. This is because the higher fre- 
quency vibrational modes (that couple to localized energetic interactions) 
are frozen out of the equilibrium distiibution by quantum statistics at low 
temperatures, and hence the rates at which the higher frequency vibrational 
modes con be excited are very small in a lattice below 1°K by detailed 
balancing considerations. 

+A hydrogen atom initially with 10 MeV energy will lose its electron and 
travel as a proton during about half the time it spends in motion in solid 
hydrogen; it is understood tacitly in our discussion that "hydrogen atom" 
means "hydrogen atom or proton" 

*^The dissociation energy of a ground-state molecule is about 4.5 eV, but 

electrons having at least 11 eV are required for the quantum mechanical 
threshold of the dissociation reaction e + H,-^-e + 2H, which involves a 

1 1 3 

transition of the molecule from the Z ground-state to a I repulsive state. 
Near threshold, each free H emerges with about 3.2 eV of translational kine- 
tic energy. 
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It Is estimated that less than .03% of the total beam power will eventually 

be absorbed by lattice vibrations*. On the other hand, the vibrational 

excitation of molecules (by impact with primary and secondary hydrogen 

atoms and free electrons) will absorb a more significant fraction of the 

total beam power, but the vibrational excitation of an Eg molecule in the 

lattice will increase the probability and lower the energy requirement 

for the chain-branching dissociation reactions e + Hj -»* e + 2H and 

H + (see ^ ru hl ar ) • It is estimated that between 25% and 50% 

ft 

of the original energy in the beam will be absorbed eventually by chain- 

branching dissociations of Hg molecules, which generate additional fast- 

moving free hydrogen atoms as well as residual free hydrogen atoms. The 

12 

latter free H atoms may become trapped inters titially (Bondybev ) or 
assume lattice sites occupied previously by Eg molecules which have been 
dissociated. If we take .375 as a tentative estimate of the energy frac- 
tion finally absorbed by chain-branching dissociation reactions, it fol- 
lows that a hydrogen atom with an initial energy of 10 MeV will instigate 
ths eventual production of about 1,66 million free hydrogen atoms In the 
solid. Conditions will be favorable for the trapping and stabilisation of 
the free hydrogen atoms If a magnetic field of about 100 toG Is applied 
through the solid. 

tlKe associated local heating throughout the solid (less than 60 1 7 * 14.3 
cal/ sec through the entire volume, giving a local heating rate of about 1.8 

XlflT^al/ cm^-aec) it controllable by virtue of the vary high thermal dlffual- 

vity of aolid hydrogen, of the order 10 ^/mc for temperatures below l°fc. 

ft 

Bnergy not absorbed by the dissociation of Bg molecules and lattice vibra- 
tions must eventually be emitted as visible and ultraviolet radiation, 
resulting from electronic recombination of ions and the radiative decay of 
vfbratlonally excited £L molecules and electronically excited V atoms and 
Hg molecules. * 



23 

In this manner about 2.0 x 10 free hydrogen atoms will be produced 
and trapped in the solid during each second of bombardment. Hence, about 
15Z of the molecules in the solid will v ‘ergo effective dissociation 
to yield trapped free hydrogen atoms during a total bombardment time of 
about 5.5 minutes. The resulting H-H 2 solid propellant wafer, 15% free 
hydrogen atoms by weight, will have a stored recoverable specific energy 
Q^c7.8l kcal/gm and thus will yield a theoretical specific impulse 
I g p & 742 sec for rocket engine conditions that make formula (A. 6) applic- 
able. 

It should be noted that the overall efficiency for energy storage in 

the H-Hj solid propellant is approximately .28, with 250 kW cyclotron 

3 

power for 5.5 minutes yielding 5.55 x 10 kcal of stored energy in a .71 kgm 
wafer of solid hydrogen. To produce a 1000 kgm amount of the propellant 
would require a conveyer belt operation and 129 hours of total bombardment 
time. The strong magnetic field must be maintained, along with the low 
temperature, for safe storage of the propellant, and the atomic hydrogen 
concentration can be monitored conveniently with electron spin resonance 
measurements during and after the bombardment of a wafer. 

* 

Appendix D: Manufacture of an H-Ha Solid Propellant 

by Impregnation with Radioactive Phosphorus 

32 

The radioactive isotope of phosphorus P decays solely by 8-ray 

32 

electron emission (to yield sulfur S ) with a half-life of 14.3 days and 

a mean energy of 690 keV for the emitted 8-ray electrons. About 300 times 

32 

greater in specific activity than tritium, the radionuclide P has been 
extensively used In medical applications and la available commercially. 
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Consider a volume of liquid hydrogen Impregnated uniformly with radio- 

32 

active phosphorus atoms In the ratio of one P atom for every 50,000 

32 

molecules (i.e., .002% molar P ) and cooled gradually through its freez- 
ing point and down to a temperature below 1°K, which is maintained for a 

32 

two-week "curing” period. During this time about half of the P atoms 
will emit 3-ray electrons having a mean energy of 690 keV and a mean range 
In the low-density solid of about 2 cm. Almost all the energy of the 
6-ray electrons will be absorbed by ionization processes which yield se- 
condary electrons with energies between 15 eV and 10 keV. The secondary 
electrons will in turn give up their energy through further ionization 
reactions, vibrational and electronic excitations of molecules, dis- 
sociations of H 2 molecules, electronic excitations of H atoms, and lat- 
tice (collective molecular) vibrational excitation reactions. At least 
a tenth of the original energy in the 6-ray electrons will be effective 
In producing free hydrogen atoms via secondary electron Impact dissocia- 
tion reactions and free hydrogen H -I- 3H dissociation reactions. Since 

4 

the net energy absorbed per dissociation is 4.5 eV, at least 1.53 x 10 
molecules will be dissociated by the secondary energetic particles that 
are generated by a 6 -ray electron with 690 keV initial energy. Thus, during 
the two-week curing period, at least 15% of the molecules will have under- 
gone dissociations to yield free hydrogen atoms. If a magnetic field of 
about 100 kG is applied to the solid, conditions will be favorable for the 
trapping and stabilization of the free hydrogen atoms in the 82 matrix. 

The resulting H - Hj solid propellant, at least 15% atomic hydrogen by 
weight, will yield a specific Impulse 742 sec for rocket engine condi- 
tions that make formula (A. 6 ) applicable. 
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Appendix E: Deflagration of an H - Ha Solid Propellant 

Consider an B - solid propellant, 15% atomic hydrogen by weight, 

at an Initial temperature of 1°K and a constant pressure of 100 psia. If 

an elemental volume of the solid is heated, the density will remain appro- 

3 

ximately constant at p * .088 gm/cm until the melting point Is reached 

6 

at a temperautre Just above 14°K. The density of the fluid, equal to about 

3 

.073 gm/cm at a pressure of 100 psia and a temperature near the melting 
point, will decrease continuously with increasing temperature according to 
the formula 

p Ss (.170 gm-°K/cm 3 ) (T — 11.67°K) _1 (E.l) 

for T > 14.0°K without significant discontinuity at the liquid-gas transi- 
tion temperature ~32°K at 100 psia (because the critical point for is 

at 33°K and 189 psia). Continuity of mass flux in the case of steady one- 
dimensional deflagration yields the relation 

P u ■ P s V (E.2) 

where u la the local fluid velocity in a coordinate frame for which the 
melting surface Is at rest and v # is the constant burning velocity of the 
solid (i.e., the rate of linear regression of the melting surface with re- 
spect to the solid at rest). It follows from (E.l) and (F.2) that 

u A .518 V s - 11.67^ (E.3) 

for a fluid element at temperature T. 

Mow with the propellant mixture yielding the overall net energy release 
Q * 7.81 kcal/gm, the functional dependence Cp« c^(T) for Kg gives a final 
(adiabatic flame) temperature of about 2090°K for the product <**• H 
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At denotes the time interval following melting for a fluid element to 
attain the temperature of 1000°K, then the temperature of a fluid element 
will Increase by the amount AT ■ 1000°K — 14°K » 986°K in travelling the 
distance At 


As = 



o 


(E.4) 


from the melting surface. Letting w denote the weight fraction of H 
atoms at any point in the fluid, we have the recombination rate equation 


du - 2 


(E.5) 


in which 


f «* 1.61 x 10 


11 


(V (u°J 


sec 


-1 


(E.6) 


is the rate function associated with a velocity-averaged total cross- 
^ ~16 2 

section of crS 1.1 x 10 cm for the recombination reaction. The tem- 
perature of a fluid element will increase from T & 14°K (for which 
u »15, see below) to T & 2090°K ( for which u ■ o) with an approximately 
linear^ dependence on u, 

« 

T A 2090°K - (13,840°K)w . (E.7) 

By putting (E.3) and dt * f“ 1 d(w“ 1 ) into (E.4) and using (E.6), (E.l) and 
(E.7), we obtain 


•Hrhe effective specific heat of the fluid is bolstered at low temperatures 
by the heat of vaporisation (~107 cal/gm) distributed about 32°X and by 
the heat of conversion of para— -Hj to ortho— Hj HL70 cal/gm converted). 
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1000 

Ax« (8.65 x 10" 8 sec) V s \ (E * 8) 

»4 

* (7.0 x 10 ^ sec) v . 

8 

Thus, the thermal gradient for heat flux to the melting surface is given 
approximately by 


~ & ( 1.41 x 10 9 °K/sec) v/ 1 


(E.9) 


The latter quantity enters the heat flux balance at the melting surface. 


X 


AT 

Ax 


v . •>. V 


CE.10) 


where X .25 mW/cm-K “ 6.0 x 10 ^ cal/cm-sec-°K Is the estimated^ thermal 

conductivity of the fluid at 14°K and h^ &■ 22 cal/gra is the enthalpy re- 

o ft 

qulred to warm and melt a unit mass of the solid Initially at 1 K . 

Hence, from (E.9) and (E.10) we obtain the burning velocity 


* 210 cm/eec. (E.ll) 

The essential correctness of the deflagration model Is confirmed by 

a 

consistency estimates. First note that an elemental volume of the solid 
will be warmed to the melting point In a time given approximately by the 


tA significant fraction of ortho-H, will be present In the solid as a by- 
product of the propellant manufacturing process, and the ortho-H- and atomlg 
H fractions will suppress the thermal conductivity from the value .fSmW/cra- K 
for liquid para-^ at 14°K and 100 psla. 

++ We have 14 cal/gm coming from the heat of fusion, 5 cal/gm from the T 8 inte- 
gration, and about 3 cal/gm from the X anomaly In heat capacity at about 
1.6 k (with a high fraction of ortho-fl, resulting in the solid as a by-product 
of the propellant manufacturing process). 
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2-9 -42 

expression D/v & 8.0 x 10 'sec, where D = 3.5 x 10 cm /sec 

is the thermal diffusivity of the solid at a temperature near the melt- 
ing point, and thus a negligible fraction of the hydrogen atoms vrill have 
time to recombine while in the solid and above the critical temperature 
(~4°K) for significant escape from trapping sites. Next observe that 
the distance Ax =£=• 1.47 x 10 ^cm [given by (E.8) and (E.ll)] has the same 
order of magnitude as the mean free path for Hj molecules at 1000°K and 

100 psia. On the other hand, the magnitude of Ax is small compared to 

-2 

the entire flame thickness, which can be estimated to be about 1.81 x 10 

cm by replacing the upper limit in the integral in (E.8) with the value 

2000. Finally note that the neglect of species and thermal diffusive 

transport in writing Eqs. (E.5) and (E.7) is justified, because the dif- 

fuslvltles in the fluid are very small compared to the mean flow velocity 

u | ^ 1.07 x 10 5 cm/sec multiplied by the approximate flame thick- 

*T-1000*K 

_2 

ness 1.81 x 10 cm. 


i 



I 
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